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ABSTRACT 


The in-line force acting on smooth circular cylinders 
placed in oscillatory flow has been measured using a recently 
eonstructed U-tunnel. 

The drag and inertia coefficients have been determined 
through the use of the Fourier analysis and found to depend 
on Reynolds number as well as Keulegan-Carpenter number for 
Reynolds numbers about greater than 25,000. A new parameter 
namely the "Frequency Parameter" has been defined as the 
ratio of the Reynolds number to the Keulegan-Carpenter number, 
and the existence of unique relations between the frequency 
parameter and the drag and inertia coefficients have been 
Shown. 

It is recommended that the experiments be extended to 
even larger Reynolds numbers uSing the same experimental 


technique and apparatus. 
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NOMENCLATURE 


Virtual amplitude of the motion 

Amplitude of the motion 

Maximum acceleration 

Average drag coefficient 

Average inertia coefficient 

Maximum in-line force coefficient 
Diameter of the test cylinder 
Instantaneous total force acting on the test cylinder 
Drag force acting on the test cylinder 
Inertia force acting on the test cylinder 
Measured force 

Gravitational acceleration 

Distance between pressure taps and mean water level 
Keulegan-Carpenter number 

Length of the test cylinder 

Reynolds number 

Distance between the pressure taps 

Pemiea Of ©Serl lations 

Time 

Maximum velocity 

Instantaneous velocity 

Width of the test section 

Frequency parameter 


Speeiuae weloghnt of fluid 








Differential pressure 
Phase angle 
Fluid kinematic viscosity 
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ie LNTRODUCTION 


Information about the time-dependent forces acting on 
bluff bodies in general and circular cylinders in particular 
has considerable practical interest in ocean and wind engi- 
neering, aS well as in basic understanding of fluid mechanics. 
Extensive discussion about the flow-induced forces and oscil- 
lations exists in the literature, but the basic hydrodynamic 
data are lacking, particularly at higher Reynolds numbers 
which are of current practical interest. 

Much of the present knowledge on the hydrodynamics of the 
oscillatory flow about the cylinders has been obtained by 
means of model tests in wave channels at Reynolds numbers 
generally two to three orders of magnitude smaller than proto- 
type Reynolds numbers. These model tests, which have been 
primarily concerned with in-line forces, have disclosed that 
the drag and inertia coefficients may be correlated with the 
amplitude of the motion relative to a characteristic dimension 
of the body and that the effect of the Reynolds number is 
obscured by that correlation. The field data have shown 
rather large scatter and no systematic correlation with either 
the relative amplitude or the Reynolds number. Understandably, 
there is considerable controversy in the literature regarding 
the effect of viscosity, free surface, etc. 

In view of the foregoing considerations, the present’ 


research program was undertaken with the following objectives: 
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(a) To furnish data, obtained under carefully controlled 
laboratory conditions, about the in-line and trans- 
verse forces acting on smooth cylinders ina 
harmonically oscillating fluid at Reynolds numbers 
up to 400,000 and 

(b) To identify the physical mechanism and parameters 
responsible for the correlation or scatter of the 
force-transfer coefficients. 

This thesis deals only with the in-line forces acting on 
smooth circular cylinders. It does not deal with ocean waves, 
non-harmonic fluid oscillations, wave and current combinations 
and its consequences, diffraction effects, free-surface and/or 
wall proximity effects, fluid elasticity or hydro-elasticity 
of flexible or flexibly supported cylinders in harmonic motion, 
or finally the interference effects between neighboring 
structural elements. Furthermore, no attempt is made to offer 
a chronological and/or critical survey of the fluid loading 
on cylinders or offshore structures. Fairly complete accounts 
in the context of wave forces are given by Wiegel (1), 

Hogban (2), and Grace (3), where an extensive list of refer- 
ences can be found. Only those works which have a direct 
bearing on the evolution and discussion of the present data 
will be reviewed, in some detail, wherever appropriate. 

In closing the introduction it should be emphasized that 
the force-transfer coefficients must be obtained experimentally 
for the analysis of the harmonic flow about cylinders may 


never be forthcoming. 
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mt. VAeeyetsyOr THE DATA AND THE FORCE COEFFICIENTS 


A completely satisfactory analysis of the resistance in 
unsteady separated flow has escaped the concentrated efforts 
of many researchers. No theoretical model can, at present, 
predict the complete force and flow characteristics of a 
harmonic motion about a circular cylinder. In the absence of 
such an analysis, the most serious difficulty lies in the 
description of the time-dependent force. Other difficulties 
arise in the description and interpretation of the history 
of the flow and the effect of the vortices and turbulence in 
the flow approaching the body at a particular time in the 
oscillation. One approximate and physically meaningful way 
around these difficulties is to assume, following Stokes' 
classic analysis of the oscillating pendulum, that the total 
time-dependent in-line force may be expressed as a sum of the 
velocity-squared dependent drag and the acceleration dependent 
inertial force, each with a suitably averaged force coefficient. 
This then is the basis of the so-called Morison's (4) equation. 


For a circular cylinder it may be written as 


= - 2a 
eee Oo Cl D plUlU + O25 7 pi Deo (1) 
For an oscillating flow represented by U = -U_Cosé pe Wwalte 
6 = 27t/T , the Fourier averages of C and C are given 


D M 
by Keulegan and Carpenter (5) as 
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21 F Cosé 
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D 4 2 
0 Pp a L D 
om 
2 U_ 4k a Siné 


2 
D Isc D 
Tl 0 Pp U. 


in which re represents the measured force, L and D_ the 
length and the diameter of the cylinder, Un the maximum 
velocity and, T the period of the flow oscillations. Evi- 
dently, Ch and Cy represent the first two terms in a 
complete series expansion of the normalized force in terms of 
the odd integers in sines and cosines. Previous studies by 
Keulegan and Carpenter (5) and Sarpkaya and Tuter (6) have 
Shown that Cy and Cu as given by equations (2) and (3) 


are the most Significant ones and sufficient to represent the 


measured force adequately. 


A. GOVERNING PARAMETERS 

The coefficients cited above will have to be correlated 
through the use of suitable parameters in order to show that 
they have some degree of universality within the range of the 
parameters encountered. A simple dimensional analysis of the 
Flow under consideration (uniform harmonic motion about a 
circular cylinder placed with axis normal to the flow) shows 


that the time dependent force coefficient may be written as: 


F 


——_—___———- = f£ (U_T/D, U_D/v, t/T) " (4) 
OMe GeULED L m = 
m 
== (kK; Re, t/T) (5) 


is 





Evidently, U_,T/D may be replaced by 2TA/D or simply by 
A/D where A represents the amplitude of the oscillations. 
There is no simple way to deal with equation (5) even for 
the most manageable time-dependent flows. The evaluation of 
the instantaneous values of the force coefficients is not at 
present feasible. Another and perhaps the only alternative 
is to eliminate time as an independent variable in equation 


(5) and consider suitable time-invariant averages of the 


force coefficients. Thus, one has 
Ch 
Cu Serco Dye Re) (6) 
CEMAX 


Even this equation (6), as simple and as idealized as it is, 
gives rise to many questions: Do the averaged force co- 
efficients really depend on both K and Re ; can one obtain 
meaningful conclusions by plotting the data for a given co- 
efficient with respect to, say, K and connecting points 
having equal Re or vice versa ; how should the experiments 
be conducted so that equation (6) yields manageable plots ; 
which of the two parameters has a more pronounced effect on 
the force coefficient under consideration ; why has there 

been considerable scatter (see, e.g., Wiegel (1)) in the field 
data when plotted with respect to either K or Re ; are 
there ranges of K and Re in which the effect of one is 
obscured by a reasonable correlation of the force coefficients 


with the other? These and similar questions have been raised 
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by many other investigators and attempts were made to establish 
Suitable correlations. The state of the art is such that the 
past conclusions and conjectures can be critically scrutinized 
only through the acquisition of data of high intrinsic quality 
obtained under controlled laborabory conditions with relatively 
simple harmonic flow situations. The purpose of such an effort 
is by no means to remove the need for actual full scale 
experiments. In fact it is to encourage full scale experiments 
and to enable those concerned to interpret and better understand 
the factors effecting the force-transfer functions. 

Past experience (5,6) has shown that the force coefficients 
obtained under controlled laboratory conditions are primarily 
functions of K at relatively small Reynolds numbers 
(Re < 25,000) and that the effect of viscosity is obscured by 
the fairly good correlation between K and the force co- 
efficients. Again, previous efforts and the reasoning based 
on dimensional analysis have shown that there is an undeniable 
effect of the Reynolds number. Thus, means have to be devised 
to delineate the effect of both K and Re or some other 
viscosity dependent parameter. 

Tt 1S apparent that the maximum velocity Le appears in 
both kK = U_T/D and Re = U_D/v . Simple rules of dimensional 
analysis state that one obtains the maximum amount of experi- 
mental control over the dimensionless variables, if the 
Original variables, that can be regulated, each occur in only 
one dimensionless product. Thus, if Ua is easily varied 


experimentally, then U_, Should occur In Only sone ones cnc 
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independent dimensionless parameters. With this hint in mind, 
Re in equation (6) is replaced by Re/K = D*/vT . This para- 
meter shall be called the "Frequency Parameter" and denoted by 


B , so that: 


= 3 (7) 


Evidently, for a series of experiments conducted with a 
cylinder of a given diameter D in water (of uniform and 
constant temperature and density) undergoing harmonic oscil- 
lations with a constant period of T , 8 is held constant, 
then the variation of a force coefficient with K may be 
plotted for constant values of 8 . Subsequently, one can 


eaSily recover the Reynolds number from 
Re = BK (8) 


and connect the points, on each 8B = constant curve repre- 
senting a given Reynolds number for a family of suitably 
selected values of the Reynolds number. Such a procedure 
SiMtnatesmene Gdithiculty of trying to draw contours of con- 
Stant K , or constant Re , or constant Cy Os Cu versus 
K or Re , or K versus Re . Suffice it to note that the 


data reported here-in shall be analyzed according to the 


relationship 


Ge 
1 


a e@ert 1 cient = iE (Ke) - (9) 
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and the Reynolds number will be used in the manner cited above. 
The power of this new procedure (new, as far as the wave force 


analysis is concerned) will become apparent later. 


aL 7, 





ITI. EXPERIMENTAL EQUIPMENT AND PROCEDURES 


A. U-SHAPED OSCILLATING-FLOW TUNNEL 

Experiments carried out by Sarpkaya and Tuter [6] in the 
initial phases of the study with small smooth cylinders at 
relatively low Reynolds numbers have proved the versatility 
and the usefulness of a U-shaped oscillating-flow apparatus. 
Thus, in an attempt to achieve larger Reynolds numbers, it 
was only natural to construct a larger U-shaped tunnel. 

Among the various designs made by Professor Sarpkaya, the 
one shown in figure 1 was finally selected for construction. 
A photograph of the completed and fully instrumented structure 
is shown in figure 2. It consists of nine modules for ease 
of construction, transportation, and mounting. Each module 
is made of 0.95 cm. (3/8 inch) aluminum plates and reinforced 
1.27 x 10 x 46 cm. (1/2 x 4 x 18 inches) aluminum flanges 
welded to the plates. The modules were assembled with the 
help of an air drying silicon rubber seal between the flanges 
of two adjacent modules and 2.54 cm. (1 inch) steel bolts 
placed 15 cm. (6 inches) apart. The inside of each module 
was precision machined so that the largest misalignment was 
about 1 mm. (0.04 inches). 

Prior to the description of its instrumentation and 
operation, a few words are necessary about the general shape 
of the tunnel. The cross-section of the two legs is 


183 x 91.5 om. (6 x 3 feet). This selection was dictated 
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Figure 2. U-Tunnel 
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Figure 3. Schematic Drawing of the 
Butterfly Valve System 
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by several considerations such as available ceiling height, 
pressures to be encountered and hence the structural and 
economic considerations, desire to obtain a virtual amplitude* 
or velocity of oscillation at least twice that of the free 
surface, period of oscillation, Reynolds number and the 
relative amplitude A/D desired, natural-damping of the 
oscillations, and the magnitude and the frequency of the 
forces. The length of the horizontal test section was chosen 
larger than twice the virtual amplitude to insure fully de- 
veloped uniform flow at the test section. Finally the two 
corners of the tunnel were carefully streamlined to prevent 
flow separation (see figure 1). This design proved to be 
more than adequate for no separation was encountered, and 
also the desired frequency and amplitude of oscillation were 
achieved. 

The auxiliary components of the tunnel consisted of 
plumbing for the filling and emptying of the tunnel, butterfly- 
valve system, and the air supply system. The plumbing 
consisted of simple piping, valves, a small pump, anda 
filter. None of these will be described further. 

The butterfly-valve system (mounted on top of one of the 
legs of the tunnel) consisted of four plates, each 4.6 cm. 


(18 inches) wide and 92 cm. (36 inches) long. A 2.54 cm. 


*The virtual amplitude is defined as the amplitude of oscil- 
lation which the cylinder-experiences. Here it is exactly 
twice the amplitude of oscillation of the free surface in 
one leg of the tunnel. 
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(1 inch) steel shaft was placed at the axis of each valve 
plate. (see figure 3) Aluminum housings supported both ends 
of the shaft with self-aligning ball bearings. A 15 cm. 
(6 inches) gear was attached to one end of each shaft which 
extended beyond the bearing. All four valve plates were then 
aligned and driven by a simple rack and pinion system. The 
rack was actuated by an air-driven piston with the help of 
a three-way valve connected to the laboratory air supply system. 
The valves, in their closed position, completely sealed 
the top of one of the legs of the tunnel. The top of the 
other leg was left open. Initially, the butterfly valves 
were closed and air was admitted to that side of the tunnel 
to create the desired differential water level between the 
two legs of the tunnel. Then the valves were opened quickly 
with the help of a pneumatically-driven three-way control valve. 
Mimordetron sect the flurd in the tunnel in oscillatory motion 
with a natural period of T = 5.272 seconds. A series of 
experiments was conducted with one of the test cylinders to 
evaluate some of the experimental characteristics of the 
tunnel. It was found that the damping of the motion is such 
that the amplitude of oscillations decreases about 3 mm. 
(0.12 inches) per cycle and about 0.15 mm. (0.06 inches) per 
cycle for amplitudes half or less than of the maximum. Thus, 
over a period of four or five complete cycles of oscillation 
at any amplitude, the amplitude, velocity, and the ac- 
celeration of the fluid changed about one percent. Evidently 


thew forced oscillations of the fluid, if such a method were to 
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be employed, cannot yield the amplitude to an accuracy better 
than one percent. Additionally in such a method one has to 
contend with some high frequency vibrations, however small, 
Superimposed on the acceleration. These result from the 
cyclic operation of the butterfly valves. It is because of 
these considerations that the experiments were carried out by 
letting the system damp out the amplitude over many cycles of 
oscillations. The advantages of the method adopted become 
apparent very quickly, Firstly, the oscillations were so 
smooth and quiet that one could not know or even hear that 
5000 gallons of water were in oscillation. The elevation, 
acceleration, and all force traces were absolutely free from 
secondary oscillations so that no Neen whatsoever were 
used between the transducer output and the recording equipment. 
Secondly, the method adopted enables one to cover all possible 
values of K for a given 8 and see the evolution of the 


forces over a period of about 30 minutes. 


B. CIRCULAR CYLINDER MODELS 

Seven circular cylinders with diameters ranging in size 
from 20.3 cm. (8 inches) to 5 cm. (2 inches) have been used 
in this study. The cylinders were turned on a lathe from 
aluminum pipes or plexiglass rods and polished to a mirror- 
Shine surface. There is no doubt that the surfaces were 
hydrodynamically smooth. The length of each cylinder was 
such that it allowed 0.08 cm. (1/32 inch) gap between the 


tunnel wall and each end of the cylinder. As will be noted 
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later, the cylinder was prevented from moving towards one or 
the other wall by means of small O-rings attached to the round 
cantilever end of the force transducers. A double-ball pre- 
cision bearing (SKF-2303-J) with an approximately 1.5 cm. 

(0.6 inches) bore was inserted at each end of the cylinder 

in aluminum housings which sealed the cylinder air tight. 

The other face of each bearing was flush with the end of the 


cylinder. 


C. FORCE MEASUREMENTS 

Two identical force transducers, one at each end of the 
cylinder, were used to measure the instantaneous in-line and 
transverse forces. The basic transducer was manufactured by 
B.L.H. Electronics, Inc. under the trade name LBP-l and 
catalogue no. 420271. A photograph of the gage assembly is 
shown in figure 4. The gage had a capacity of 2224 N 
(500 lbf.) with an overload capacity of 200 percent. The de- 
flection of the gage under 500 1lbf. load was 0.25 mm. (0.01 
inch). For the largest cylinder and amplitude encountered 
in the experiments, the maximum load was less than 200 l1bf., 
and the deflection of the beam was less than 0.2 mn. 

(0.008 inches). 

A special housing was built for each gage so that it can 
be mounted on the tunnel window and rotated to measure either 
the in-line or the transverse force alone. The bellows which 
protected the strain gages had to be waterproofed in such a 


manner that they would not adversely effect the operation 
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Gage Assembly 


A. 


Figure 





of the gages when subjected to about 6 m. (20 feet) water 
PeecoouUre ae temperatures 18 C (64 FP) to 74 € (165 F). For 
this purpose the bellows were completely filled with Dow 
Soring 3140 (RTV coating without bringing the rubber into 
contact with air during the filling operation. After filling, 
the ends of the bellows were sealed air tight with special 
clamps. The silicon rubber remained in its original liquid 
form throughout the operation of the gages. 

The cylinders were placed in the test section by re- 
tracting the gages from their housing and then pushing them 
into the bearings mounted on each end of the cylinders. As 
noted earlier the O-rings placed on the cantilever end of each 
gage prevented the test cylinders from moving sideways to- 
wards one or the other wall and helped to exactly set the 
0.8 mm. (1/32 inch) space between the cylinders and the 
tunnel wall. The cylinders were free to rotate at the 
Bpplreation Of a slight torque by hand. 

After the mounting of the first cylinder, the exact 
angular position of the gages within their housing had to 
be determined and set with a pin so that the gages measure 
either only the in-line or the transverse force. For this 
purpose, an approximately 400 N load was hung on the cylinder 
with a lubricated nylon rope. The in-line force (acting in 
the horizontal direction) was observed on the amplifier re- 
corder system. Then the gage was rotated in small increments 
until the in-line force was exactly zero. A final check was 


made by measuring the outputs of the gages with a precision 
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voltmeter. The position of each gage was marked and set with 
a pin. Finally, four bolts were placed on the gage housing 
to hold the gages rigidly in position. Removal of these 
bolts and the pin allowed the rotation of the gages exactly 
90 degrees, after which the bolts and pin were placed in 
position. In this manner the gages were capable of measuring 
either the in-line or the transverse force without any "cross 
talk" between the in-line and transverse forces. Ordinarily, 
one gage was set to measure in-line force and the other gage 
the transverse force. At times both gages were used to measure 
only the in-line or the transverse force. 

The calibration of each gage was accomplished by hanging 
loads in the middle of the cylinder after setting both gages 
to sense only the transverse (here vertical) force. The di- 
rectional sensititivity of the gages was also checked by 
applying identical loads upwards on the test cylinders with 
the help of a hook-cantilever arm attached to the top of the 
tunnel outside the test section. Repeated calibrations have 
shown that the gages are perfectly linear up to 2000 N; they 
yeild the same signal for loads applied either downward or 
upwards; and that the gages, together with the electronic 
system to which they were eventually connected were capable 
of sensing forces as small as 0.1 N (0.02 l1bf.). 

The analog data appearing on a two channel Honeywell 
recorder in the form of either in-line force versus amplitude 
or transverse force versus amplitude were were read at every 


0.1 seconds. The evaluation of the data is discussed later. 
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D. ACCELERATION, ELEVATION, OR VELOCITY MEASUREMENTS 

It 1S because of the extreme importance of the accurate 
measurement of the instantaneous value of these quantities 
that they are discussed here separately. 

Firstly, it should be noted that the measurement of the 
amplitude, acceleration, elevation, or the velocity is a 
matter of interpretation of the Signal received from the 
appropriate transducer in light of one of the following 


expressions 


tw ieeh T 5.272 seconds for the experiments reported here. 
Three transducers were used to generate three independent 
d.c. Signals, each proportional to the instantaneous value 
of one of the quantities cited above. The first one con- 
sisted of a seven foot long platinum wire stretched vertically 
in one leg of the tunnel. The output of the capacitance-wire 
bridge was connected to an eight channel amplifier-recorder 
system. The response of the wire was found to be perfectly 
linear within the range of oscillations encountered. The wire 
was capable of yielding a meaSurable signal for changes in 
water elevation as small as 0.8 mm. (1/32 inches). Such a 
sensitivity was not, however, always desirable for the in- 
Stabilities on the water surface gave rise to small oscillations 


in the analog records. The effect of such instabilities were 


practically eliminated by placing the wire along the axis of 


Zo 





SOc. fLOOtL) sdiameter and 213 cm. (7 feet) long thin 
plastic pipe. Be that as it may, the use of this wire was 
rendered unnecessary due to the use of a more reliable method. 
The second method consisted of the measurement of the 
instantaneous acceleration by means of a differential-pressure 
transducer connected to two pressure taps placed horizontally 
bl cm. (2 feet) apart and 91.5 cm. (3 feet) to one side of 
the test section. The output of the transducer was again 
connected to the eight channel recorder. The instantaneous 
acceleration was then calculated from Ap = ps(dU/dt) where 
Ap is the differential pressure, s the distance between the 
pressure taps and dU/dt is the instantaneous acceleration 
of the fluid. The effect of the pressure drop due to the 
viscous forces over the distance s was calculated to be 
negligible. 
The third method again consisted of the measurement of 
the differential pressure between two pressure taps. The two 
taps were placed symmetrically on the two vertical legs of 
the tunnel at an elevation 96.52 cm. (38 inches), (see 


figure 5) below the mean water level, i.e., H = 38 in. 


re instantaneous level 


A, 
mean level am ic 
Ny 





pressure taps 


Figure 5 
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Applying Bernoulli's equation for unsteady flow between each 
pressure tap and the instantaneous level of water, it is easy 
to show that twice the amplitude of the free surface oscillation 


(virtual amplitude) is given by 


i = Ap/y 
l Sey ale H 
oom k 


in which g and T are constant and H is kept constant. 

Thus the signal of this transducer yielded the virtual amplitude 
or the maximum velocity in each cycle. It was entirely free 
from noise or small free surface effect. The transducer was 
calibrated and its linearity checked before each series of 
experiments. 

Suffice it to say that all three methods gave nearly 
identical results and yielded the amplitude, velocity, or 
acceleration to an accuracy of about two percent relative to 
each other. These comparisons as well as the perfectly 
Sinusoidal and noise-free character of all pressure and force 
traces speak for the suitability of the unique test facility 


Weed in thas study. 


Ba Data REDUCTION 

Experiments were repeated at least three times for each 
cylinder. Often only two, but at times as many as four runs 
were evaluated. Suitably selected cycles were read every 0.1 
second and the tabulated data were transferred to the computer 


data cards together with additional relevant information such 


ot 





as calibration factors, cylinder diameters, number of data 
cards, etc. A computer program calculated the force transfer 
coefficients previously discussed (see Appendix A). 

Sample traces of the in-line force and the elevation are 


shown in figure 6. 
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Elevation and Drag Force Traces 


Figure 6. 
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i ee SCUSSION OF RESULTS 


A. BLOCKAGE AND LENGTH-TO-DIAMETER EFFECTS 

Attempts to achieve as high Reynolds numbers as possible 
in conducting wind-tunnel or water-tunnel experiments in- 
variably give rise to wall-interference effects which, of 
course, influence whatever measurements are desired. There 
are several blockage correction formulas for steady flows 
which might be used so that the wall-interference effects 
might be minimized. Unfortunately none of these formulas 
could be used in the present study for no one has demon- 
strated that the blockage effects in oscillatory flow are 
identical to those experienced in steady flows. 

The blockage ratio D/w , where D is the diameter of 
the cylinder and w the width or height of the test section, 
and the length- or span-to-diameter ratio, L/D , for the 
cylinders used in the present study are tabulated below 


oe — 99 .44 cm. (3 feet), L = 90.885 cm. (2.9818 feet)). 


D D/w L/D D*/vr 


20.244 cm. (7.970 inches) Onze 4.49 1908.7 

16.447 cm. (6.475 inches) 0. 5.52 525929 

Bei 7ecme (52975 inches) O'2 Ag Bb 9 4480.2 
12.674 cm. (4.990 inches) 0.14 7.17 3123.2 

10.103 cm. (3.978 inches) Oasalgal oe BS, 1985. 
7e544 ,eme12.9/70 inches) On0e2 ie 20 5 BUG. 
6.349 ‘om. (2.4996 inches) 0.069 Panes ySae 

0 


5.057 cme. 991 1nches) s055 IS ADF. 


No co dD NH 


5 


* y= 2.05 x 10> ft*/sec (0.009755 stokes) 
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Also shown in the above table are the 8 = D?/vt values for 
later use. Achenbach (7) and in some of the experiments of 
Fage and Falkner (8) the blockage ratios were 0.166 and 0.185 
respectively. The length-to-diameter ratio in Fage and 
Warsap's (9) experiments was 20.2 or 7.88, depending on the 
diameter of the two cylinders they used, as compared to 3.33 
in the experiments of Achenbach. It is generally observed 
that values of Ch are smaller for cylinders with larger 
length-to-diameter ratio. Thus the presence of both the gaps 
amas the larger L/D ratios could result in lower drag 
coefficients. Furthermore, the wake of the cylinder, near 
the cylinder ends is supplied with high pressure fluid from 
the front and as a result smaller values of Ch are expected 
Since the base pressure iS increased over the value it would 
otherwise attain. In the subcritical range of Re, however, 
these effects appear to be negligible. In fact, Mosbach (10) 
found that in the subcritical range there is no effect of 
length-to-diameter ratio. 

In the present experiments, the gap cannot be eliminated 
by extending the cylinder into a cylindrical cavity within the 
two windows supporting the gages and the cylinder because of 
the fact that eight cylinders of different diameters were used. 
Obviously, it would have been too costly to build a window 
for each cylinder. It is believed that the very small gap 
used plus the cantilever end of the gage extending into the 
bearing in the cylinder minimized the supply of high pressure 
fluid into the wake from the front of the cylinder during one 


half of the cycle. 
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Returning to the discussion of the blockage effect it 
must be emphasized that the formulas used for steady flow 
correction effects cannot be applied to oscillating flows 
and that there is not a unique blockage correction for the 
entire period of the harmonic flow. This is evident from 
the fact that within a given cycle the fluid undergoes varying 
accelerations and velocities as the wake width, momentum 
deficiency, and the wake pressure change accordingly. fThus, 
a blockage correction made for the instant of maximum 
velocity is not applicable to the instant at which the maximum 
acceleration occurs. In view of the fact that there are no 
previous investigations, a series of experiments had to be 
conducted to determine the role of blockage in the flow under 
consideration. For this purpose a differential pressure 
transducer was connected to two pressure taps on the same 
Side of the tunnel wall. One of the taps was placed on the 
wall directly above the axis of the test cylinder. The other 
tap was placed 76 cm. (29.92 inches) to one side of the first 
tap along a line parallel to the flow. A Series of experiments 
were carried out with the 16.447 cm. (6.475 inches) cylinder 
and the differential pressure was recorded and compared with 
the differential pressure obtained from the acceleration 
transducer. Furthermore, to simplify the comparison both 
transducers were calibrated so as to render exactly the same 
output under identical calibration loads. The results have 
shown that the two differential pressures were nearly identical 


and that they were certainly within three percent of each other. 
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Often the two traces of two transducers were undistinguishably 
coincident. This somewhat surprising result is a clear 
indication of the fact that the blockage effect in harmonic 
flows 1s negligible at least for D/w ratios less than 0.18. 
Although no special attempt was made to interpret the lack 

of blockage effect in such flows it is believed that the 
presence of vortices on both sides of the cylinder together 
with the high periods of acceleration and velocity render 

the flow relatively more uniform at short distances away from 
the cylinder in the test section. Therefore, for the reasons 
cited above no blockage-effect corrections were applied to 
the data presented here. It might be of interest to note 
that had the flow been assumed uniform and had the maximum 
velocity for the largest cylinder and the Reynolds number 
were used to calculate a blockage effect correction through 
the use of one of the exiSting formulas, one would have 

found that such a correction would have amounted to about 


ee percent. 


B. IN-LINE FORCE DATA 
In-line force coefficients calculated as previously de- 


scribed are plotted as a function of the so-called Keulegan- 


Carpenter number denoted here by K . A representative plot 
of Cy wekhsus K for the 16.447 cme (62475 inches) cy Jandex 
HomSheowein figure /. It 1s evident fromythis figure that 


there is practically no scatter in the data even though it 


represents the result of four independent runs. This is 


37 





once again indication of the fact that the tunnel and its 
instrumentation are extremely reliable and that the data when 
evaluated with extreme regard relative to the phase angle are 
perfectly repeatable. As noted earlier the Reynolds number 
along this plot varies with K and 8. One can designate 
special values of the Reynolds number along this curve simply 
by writing Re = Kf§ or by finding the K values which 
correspond to a particular value of Re by writing K = Re/8 
Prior to doing so mean lines have been drawn through the data 
as carefully as possible for each cylinder data. Then round 
numbers such as 30,000, 40,000, etc. were chosen for Re and 
the corresponding K values were calculated and indicated on 
each plot. This simple procedure resulted in figure 8 in 
wach C is plotted as a function of kK for all cylinders 


D 
and the identical Reynolds number points are suitably 


connected. This procedure has also been repeated for the 
inertia coefficient Cu and the result is presented in 
figure 9. It is evident from figures 8 and 9 that there is 


a remarkable correlation between the force coefficients, 
Keulegan-Carpenter number and the Reynolds number. The 
smoothness of the constant Reynolds number lines is another 
indication of the consistency of the data from one cylinder 
to another. 

Figure 8 shows that for a given Reynolds number, for 
Reynolds numbers larger than approximately 50,000, Ch 


increases at first with increasing K , reaches a maximum, 


then decreases and then continues to rise gradually. For 
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Reynolds numbers larger tnan about 50,000 thesinitial rise 


lige - with K is very small or non existent. For very 


large Reynolds numbers and K , Ch Aten S scmpmwacilecdi hy, 


constant value of about 0.5. For a given K , however, Cy 


always decreases with increasing Reynolds number. It is 


evident from figure 8 that the change in Ch with Reynolds 


number for Re smaller than approximately 25,000 is almost 
negligible. It is for this reason that Keulegan and Carpenter 
(5) and Sarpkaya and Tuter (6) reached the conclusion that 


Cy does not depend on Reynolds number. The variations of 


Cp with Re become most apparent between Re = 25,000 and 


Re = 150,000. It is also noted that for large values of 8 


Such as 8 = 5260 and 4480 the Ch values do not significantly 


change with Re . In other words Cy values corresponding 


to 8B > 5000 yield the lower limit of the Cy values whereas 


Cy values corresponding to 8 = 500 consitute the upper 


mimet, of the Ch values. The most significant changes in 


Cy with K and Re occur between these two 8 values. 


Figure 9 shows the variation of the inertia coefficient 


Cu with K and Re. The most significant changes in C 


occur in the neighborhood of K =12. It is also around 


M 


this value of K that C changes Significantly with Re. 


M 
As discussed in connection with the variation of the drag 
coefficient, 8 values of 500 and 5000 more or less de- 
termine the lower and the upper limits of the inertia 


coefficient. Figure 9 also shows that Cy for a given K 


increases with increasing Reynolds number and varies very 
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little for Re smaller than approximately 25,000. It is 

for this reason that in the studies by Keulegan and Carpenter 
(5) and Sarpkaya and Tuter (6) Cy was thoughtto be inde- 
pendent of Reynolds number. The present results show that 
this 1s not so and that Cu varies with both Re and K 
throughout the K and Re range encountered in the present 


investigation. It appears that for very large Reynolds 


numbers Cu will approach a value of about 1.85. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The present investigation of the in-line forces acting on 


smooth cylinders in the range of K values from about 3 to 


200 and Re from 10,000 to 400,000 warranted the following 


conclusions. 


is 


The drag and inertia coefficients could be determined 
quite accurately through the use of the present 
experimental system. 

Both the drag and inertia coefficients depend on the 
Reynolds and Keulegan-Carpenter numbers. For a given 
K , the drag coefficient decreases with increasing 
Reynolds number whereas the inertia coefficient 
increases with increasing Reynolds number. 

The most significant changes in both Ch and Cu 
occur in the neighborhood of K= 12. 8 vaiues of 
about 500 and 5000 determine the upper and lower 
limits of the two coefficients. 

Neither Cy nor Cu vary appreciably with Re for 
Rew 825,000. 

It is recommended that the experiments be carried out 
at even larger Reynolds numbers through the use of 
the existing systems by decreasing the viscosity of 


the Water through heating and with cylinders of 


various artificial roughnesses. 
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It is also recommended that additional investigations 
such as the hydroelastic oscillations of flexibly 
mounted cylinders in harmonic flow and the hydro- 
elastic oscillations of ridigly mounted flexible 
cylinders be carried out. Finally the determination 
of the effect of wall proximity on the in-line and 
transverse forces acting on cylinders and other bluff 


bodies will be of extreme importance. 


45 











et cere 


ON eee) NaS AOAC OAT 


ane 


Cine 





— owen wee ee 


oy 


Fee) 
eal 


cs ee ee _~ 


PVA DON pcr 4 


COC YO Ga ee aes 


Maan | Te 0A te TGs oan 


APPENDIX A 


COMPUTER PROGRAMS 
re cee COCEe CECE CC CCE CCPC CCE CCCE CC CC CCCC CEC CCC CECE 


sewn Ay or, 


— «= NR i a er a cert ee 


MVEEACE CD Al DCm CALCULA TSCNIONEOR CYL i” 


ar: ie 
(UMAX*PLR/DI A) 


ee =O ye ihe Ss 3 ae 
CER =ODEMCMS CE NLE ES MISPL 
Me eer Li ue a OM UTR 
poe t+ CoEPrit Pe. Ss 
eehee eA CUOEFE ICI: hm 
tP= PER FID Tie SECONDS (CHART 
Catv? i> ier PO Eh 
Sf ae Ae OEE ob UNC aN 
Toe rete “CSL FY Urea 
we iclivwe. ot @e. SES 
We Cee? Deh) Beer &, 
Pee (ae Gl CCC CL CE 
rll FOR CE{( 36) 
iP LOD. ah 
=e 6 a 4 : 
Pr Soy A pe CK ys Keds SG 
'=3,14159 
Poi tes. Few 5 
L=1.4946 - 
ee (5, EA et Py fee, cay a 
Pa eo tee) / (Pye ps OF 
LZ =(=3=PEASPEP / (OXF Hi mL? 
Ie eA 
Coe Pur tr Be A 
Dee 2 ‘ 
tee et cape eA) / (MU 
WOME (6,15 F011 A ee fs 
mre ECG. Zi) 
lL =C, ) 
ere (et 
aay & KeL POS RD 
aU) 


a 
as 


<a P| 
AAG 


= 


tes 5 


ale We aaal 


t eeTI oes ALN: 


Sos (C2 haa 
——— 
bee 


<= 
See 
~~ 
« 
> woe 
. 


> one 


em 
te 


a Maen CC Cae CS Nath 


TY te ee 


cq 
e868 SCODO 


TN Cr ie Ss. Cee GI 
"OG Do 


MMS COC arr tt it 


© 


= 
OV 


FL Sey EAS em NN SRE SOLAR! pe 


-C, 

=. 
eee) O27 91202 
eee eek =i, §.C4+D 
=C{ *FCFCECK ) 
Ae 2 Ps ae 
Mies 1 (| Oidh } 
wot —-Ce. (Al PAL) 
Sli [at TAT gm te 
Clea we FSA 
Petite rat Lear SLCl 
Peo eee Pee DEL eT Sie as 
BBSHE22PplxFeruc SA*ABS(C 


Dl 
Yi D> 
Ph 


nee CCCECLECE CLE CELL. 


Re tee ene Me he et Acie - 


ee 1X CN, aan 


aes ea 


rage eet vie 
rims co Uae fol eC Lx AUPE AMD) 


~ 
wh 


—2 


PER-IGD/ChA~—T— SPEED) 


COCCECCECUCS CC 


Oe ee el 


MP, PER ,UMX, CN NE UNE 


| 
: 
i 
| 
| 


CVA VO MIEN ICC? =O Gea 


Seana 








{ | 
. rm 
| O 3 | 
| O | | | 
: | n ~*~ 
, Ye { ! 
Le 1 
| W | 
H yee ' 
] | os 
<{ . 
a) | i) 
| Cai. | | uw 
Lit 6 Es | us O 
| Eee | | QO. 
<I 3 =) SS _- 
27 Oe } ~- 
| | a Un : | a ! ~ Hl 
| | oe) Lc | | C% a Pe 
: | Le ti oS ~ 
os ° + 2 : my . i ae 
it ' «a “Ti | | mae 
oe iS) Li ~ | — | iv | Ge 
2 | C3) 7 ae = v2 | .. 
b | | oO cm Ff : 3 ul | on 
i | ee co | C | va | -N 
idJ | - ok, |e t A o- ioe 
es : ~ > a a n-t  T ss aa 8 
— | <i j mm) WJ =, | | Os <1 0 
<z UL | i [eet a toe | 4 ert: 
Ur 3 : a) Mize OO Poe vy ~o 
ey; | | oa) a ~ Tts as > 2 | ar 
C) <i ~ oe Se «fifa a mes <TLL : Qu 
wena. Se <Z 0 me aad tend “Mae ¢ o, 
ce Bt iw. a oad oh ae ~~ ‘eo We ona 2) 
LN « | a] Lal am aa, ae i o= 
I Chk «> [—~ ~~ es eee Loi | ae 0s 
3¢ 3h - io ~ wit Ge = “on Of WU * 
<“T <1 | | | < <\ ws <j —_~——™ e/— ~~ a~ fl, i ol OS 
Jj. OC | pee, ae = | CoS OI ort mu rary o- © 
Cy met 9 | | o. a a | PS ee ou oa out 
Lesh 2) iC! 6 ia CNet tL ot wilh «fh u.<i e 
St ode 2 : ef lag: “© ~ “951 Te) de bie MU. SO Oya 
‘oO : ) ~ ea ee (S =) 8:-2. = a Laer a 
Pe ee — Iie Se Oe) Ge. werlloe.j~ tae, (AIC-)) TAO LE 0 <{ Deo 
bearer St x - | -- buyin oe | Go Uhwast) pa) Oa oo [ Ub bealee a oak any ri «2 
CR ek tee 2) oe) 5, bo 4. <- Te (oS NS ff a4 o> i 
ye, ee ea Saas ‘iia |S ae XSI 2 i ert iaed | a 1 m—ye + WI af OOS 
Ee a tT CL) ele yo Cae) ge ee ee ee LE mL) an UL me 1 ee... 
Hee Oe 8 NTUNO LN | AUICKS DANY ek Udi NL + ‘o= = 
reer eth bh et + a hd +. + Ott! Om™ Fat 1 OD eC) meri C7 LL , oe ee Lite 
Cid ot Ma ts aye. om (7) ele” <8 ete ee 8 Ge Hears Oo. oD itay «LS Lins 
oe at <i) Qua Ri JIS DOG HAIN ImQOSO FF Cl D4 >- at fey Be ie Sek “ue —  — o< 
QIN OUT OO TL UL eB NN Ped ee mee ee ee 4c oe eo EN tai, meme hb Tork bk ON 
woyou LiL Ln UP LS ON th Lydd i 2 4k Ph, irs ta WCW MH UN Ll eee a 
Crime?) Veet) eee LL, ne, ed ee ee, ees If LL I CI) ee 
OOiL< woud Li lf leat Ft) il pa I ee et i) oll a Bh eed BOY ee Co Ee 


ye Fey eee ye: Uk Kany a Com eC. ~ OFF O23 Otero.) ti THN et LL. Woe Usk Ge ae C LC} 
Se ei EE LD Srey. 4 Bg Gyo Ce at Swe OW WL bu eiiitco Bt het Odkk tL 
W { mi 

(a) { j CPaae at 

| @) OOD mt 
aN) 3 Ve 


| Veet vee = 


~ 

j »< 

D 
«f ~ ~ 
Pa! x< he 
oot ig Yi 
Css ~ = 
- - re 
-~ ff - 
~< C.) o 
uN a. * 
— > on) 
~ EST ~ 
- iY ~ 
A 2 i 
> ~ = 
Cy >< ‘ 
C9 i Go 
-= os = 
~< wu pe 
od Su. ©: 
- we, CD = 
oe a= 86 
ae = =) i 
G >< 
a oe - 
<Ie fh KK 
=~ -[ *® ee en 
om i Niee fi ~——_— 
Mee emt -— = 
acyl a ee cere 
|) aie. 2m? OOo 


ee we Ot Ht 

cy ~f w= ex: e 
LL ed etd De UV 
QUO R <a ew eet 
Vie we om) ON 
PG | Pio Mr 

em Uyst ts i fer at 

| tLe tig @, - e 

- oe J Cle ON 
ee ye ott 
ep oe em i EE tS 
~ 2C - &:& & & & & ~ 6 
= f= = - <e - = ee! 
SD --Or OCxXO-! 
oF = ne & “\J- te. 
set | ee et ee eet ge | eet et 
roe ie iene 
< edte1Ly)<(<fe hi Shae 
| he 3 aS s)\ 3 * oo SS. 
CONG a OG Ch Wits “ace 
WO Cs aC tlr- & 
Wwe thie Lib Wil 
ae” ales CNN LUN 

ASE OOO ste WU 


| 


—— 


' 
} 


eS 


47 





ioe 


Est, OF She BeRENGES 


Wiegel, R. L., Oceanographical Engineering, Prentice Hall, 
HG. ,. ENnGlewood =elittis, N.J.» 1964. 


Hogben, N., "Fluid Loading on Offshore Structures, A 
state of Art Appraisel: Wave Loads," Maritime 
Technology Monograph, National Physical Laboratory, 
England, 1974. 


Grace, R. A., "Wave Forces on Submerged Objects," 
Miscellaneous Report No. 10, Univ. of Hawaii Look 
Cab-M-—10,, July 1974. 


Merrson, J. R., O'Brien, M. Pos) Johnson, Joalw., and 
Schaaf, S. A., "The Force Exerted by Surface Waves on 
Piles," Petroleum Transactions, A.I.M.E., Vol. 189, 
po. Me 9—- tb 547 950. 


heulegan, G. 48. and Carpenter, Umi Forces On Gy rnecrs 
and Plates in an Oscillating Fluid,” Journal of Research 
of the National Bureau of Standards, Vol. 60, No. 5, 
Research Paper No. 2857, pp. 425—44@pemay 1958. 


sarpkaya, T. and Tuter, O., "Periodic Flow About Bluff 
Bodies, Part 1: Forces on Cylinders and Spheres in 
op oinvcotmeativ Oscillatmngahiuia.. “NPS= Technical 
Report No. NPS~59SL74091, September 1974. 


Achenbach, E., "Influence of Surface Roughness on the 
Cross-flow Around a Circular Cylinder," Jour. of 
Piven Mechs., Vol we46,. pp.) 321—-s35,-197 1. 


Fage, A. and Falkner, V. M., "Further Experiments on the 
Eltow Around a Circular €yibinder? "Acro. kes. Comm. 
London, Reports and Memoranda No. 1369, 1931. 


Fage, A. and Warsap, J. H., “The Effects of Turbulence 
and Surface Roughness on the Drag of a Circular 
Cylinder," Aero. Res. Comm., London, Reports and 
Memoranda No. 1283, 1929. 


Morsbach, M., "Uber die Bedingungen fur eine Wirbelstrassen- 


bildung hinter Kreiszylindern," Dissertation, 
iene Aachen, 1967. 


48 





0, 


INETITAD DISTRISULION Liss 


Defense Documentation Center 
Cameron Station 
molexandria, Virginia 22314 


Pioeary ,. code U212 
Naval Postgraduate School 
Monterey, California 93940 


Department Chairman, Code 59 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


Pae@ressOn sl. Sawpkaya, Code soo. 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


Lieutenant Savas Onur, Turkish Navy 

havin milo KK, Makina UZ. VE we i se@.G.e Onatan 
Keeaecli, Golcuk 

TURKEY 


Istanbul Teknik Universitesi 
Makina Fakultesi 

Taskisla, Istanbul 

TURKEY 


Dz, Hrb. Ok. Kutuphenesi 
Heybeliada, Istanbul 
TURKEY 


Ortadogu Teknik Universitesi 
Makina Fakultesi, Ankara 
TURKEY 


Neil Collins 
Naval Postgraduate School, SMC 2655 
Monterey, California 93940 


Olver 7 Ab bhesig 

Dz, Hrb. Ok. Ogretim Uyesi 
Heybeliada, Istanbul 
TURKEY 


49 


No. 


Copies 


Z 


























Onur 
In-line forces 
actina on smooth cvl- 
inders in harmonic 
f low. 


